Dark state cooling of a trapped ion using microwave coupling 
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We propose a new dark-state cooling method of trapped ion systems in the Lamb-Dicke limit. 

With application of microwave dressing the ion, we can obtain two electromagnetically induced 
transparency structures. The heating effects caused by the carrier and the blue sideband transition 
vanish due to the EIT effects and the final mean phonon numbers can be much less than the recoil 
limit. Our scheme is robust to fluctuations of microwave power and laser intensities which provides a 
broad cooling bandwidth to cool motional modes of a linear ion chain. Moreover, it is more suitable 
to cool four-level ions on a large-scale ion chip. 

PACS numbers: 37.10.Rs, 42.50.-p 


I. INTRODUCTION 

Cooling a single trapped ion to the ground state has 
been a fundamental subject in quantum metrology, quan¬ 
tum simulation and quantum computation [l|-|4|]. Side¬ 
band cooling is the first method to reach the vibrational 
ground state of a single trapped-ion and has been imple¬ 
mented experimentally |5H7|. Two requirements should 
be satisfied for sideband cooling: ( 1 ) the ion should be 
localized to dimensions much smaller than the optical 
wavelength (Lambe-Dicke limit, LDL), and (2) the ab¬ 
sorption sidebands should be well-resolved (strong con¬ 
finement) ,i.e., the vibrational frequency should be much 
larger than the natural linewidth T of the optical tran¬ 
sition for laser cooling. Although the sideband cooling 
can cool the ion close to the ground state, the final mean 
phonon number is still higher than the recoil limit due to 
the off-resonant carrier and blue sideband transitions. 

To prevent the off-resonant carrier transition, quantum 
coherent effects have been exploited in several schemes, 
e.g., the cavity-assisted cooling scheme [§0 , and the 
dark-state cooling schemes [lll4l8j . The electromagneti¬ 
cally induced transparency (EIT) cooling is the first dark 
state cooling proposed by G.Morigi et al 0 which has 
been a standard experimental scheme to cool trapped 
ions 0 , 0 ] , and cooling a nanomechanical resonator to 
its ground state has been proposed based on EIT cooling 
IS mi- EIT cooling uses a three-level trapped-ion sys¬ 
tem coupled by two lasers with the same blue detuning, 
a Fano-like absorption spectrum with a zero at the car¬ 
rier transition and a peak at the first-order red sideband 
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transition. Therefore the EIT cooling scheme can achieve 
a lower temperature than sideband cooling. Moreover no 
strong confinement is required due to the three level na¬ 
ture of the system. 

Further improvements in approaches to reach the 
ground state are achieved with the double-dark state con¬ 
cepts by eliminating all heating transitions. Accordin g to 
the EIT cooling scheme 0] , the double-EIT cooling Qj} 
results in a double-dark state for the two EIT structures, 
thus eliminating both the carrier and the blue sideband 
transition. An alternative approach that combines the 
methods of EIT and Stark shift cooling [0 is the ro¬ 
bust and fast cooling scheme [ 0 , which cancels both 
off-resonant excitations due to the combination of EIT 
and Stark shift Hamiltonian contributions. Dark-state 
cooling using standing waves 0 localizes the ion at the 
node of the laser to eliminate the carrier transition, and 
the blue sideband heating is prevented by the EIT effect 
[0 . Ground-state cooling by quantum interference path¬ 
ways 0 and fast cooling through magnetic gradients 
[181 ] can also cancel all the dominating heating effects. In 
the above schemes, all the leading-order heating effects 
are suppressed, and the whole system can be cooled to a 
steady state, namely, the double dark state. 

|4-) = |a 0 )|0)+/3|a 1 )|l) + O(r ? 2 ). (1) 

where |oo), |ai) are the atomic internal states, 10 ), | 1 ) 
are the vibrational ground state and lowest excitation 
respectively, and /3 is a coefficient of the order 0(r)). 
Thus, the final phonon number is zero in the leading 
order of the expansion in the Lamb-Dicke parameter 77 , 
i.e, 

n ss = 0 + 0 (r f) . ( 2 ) 

The scheme presented here also eliminates the carrier 
and blue sideband terms completely (up to the first-order 




in the Lambe-Dicke parameter), using an atomic four- 
level system, a combination of two laser beams and a 
microwave field. The application of a microwave field is 
used to dress two bare atomic ground states. By tuning 
the Rabi frequency of the microwave field to the condition 
fl = —1//2, two EIT structures show up to cancel both 
the heating transitions. 

Our scheme has several advantages in experimental im¬ 
plementation: (1) similar to the double-EIT cooling 10], 
EIT cooling [III, [j_(j and sympathetic EIT cooling 19], 
our cooling scheme requires relatively small laser intensi¬ 
ties and provides a broad cooling bandwidth such that 
multiple motional modes of a linear ion chain can be 
cooled simultaneously compared to previous implemen¬ 
tations of conventional Raman sideband cooling. How¬ 
ever, the scheme here uses a microwave instead of one of 
three laser beams in the double-EIT cooling, and mi¬ 
crowave manipulation components are inherently sim¬ 
pler, more stable and scalable compared to lasers [2J-|28|; 
(2) for many cooling setups [3 H3, HE HH based on a A- 
configuration, they requires additional repumping field- 
s because the excitation state decays to more than two 
lower states while the repumping process can disturb cor¬ 
responding double-dark states, and worsen cooling per¬ 
formance. Instead of repumping the system arbitrarily, 
we use a suitable microwave to improve cooling efficien¬ 
cy. Therefore with our scheme the unwanted additional 
decay pathways can become an advantage. In addition, 
unwanted upper levels can also affect the cooling result 
seriously [la, lla] , while it has a limited affect to cool the 
atoms with nuclear spins by using our cooling scheme. 

The rest of this paper is organized as follows. In Sec. I- 
I, the cooling scheme and corresponding Hamiltonian are 
presented. Under the appropriately chosen microwave 
Rabi frequency and lasers’ detuning, the system can stay 
in a double-dark steady state and all the main heating 
effects are cancelled. In Sec. Ill, the analytical formula 
for the cooling rate is obtained by the method of adiabat¬ 
ic elimination and the optimum condition for the cooling 
process is derived. Sec. IV provides numerical simula¬ 
tions of the cooling dynamics and the robustness against 
the laser and microwave intensity fluctuation. Finally, a 
brief conclusion is given. 


II. COOLING SCHEME AND HAMILTONIAN 



FIG. 1: (Color online) The four-level atom has an excited 
state |e), and three ground states | g), |r) and |d). |e) | g) 
and |e) |r) are coupled by the two laser beams while | g) o 

|d) is driven by a microwave field. 

and f l g , respectively. States | d) and | g) are coupled to a 
resonant microwave field with frequency ujmw and Rabi 
frequency Hmw■ The Hamiltonian of the whole system 
reads (H = 0) 

H = vtfb - A g \e) (e| - (A g - A r ) | r) (r\ 

+ fig |e) (g\ e lkL ^ x + fl r |e) (r| e lkL ^ x (3) 
+ Qmw | g) {d\ + H.c, 

where A g , A r are the detunings of the two lasers, respec¬ 
tively , which are 

A g Cl^L\ ^g) , 

A r ^L-2 av), (4) 

A gr = A g — A r . 

The cooling process occurs in the Lamb-Dicke (LD) limit. 
Thus, the Hamiltonian ([3]) can be expanded by the LD 
parameters around the trap center: 

HlD = Hat. + Hm + V, (5) 

where 

H m = vtfb (6) 

is the free Hamiltonian for vibrational degrees of freedom 
(DOF), 


Our scheme is shown in Fig. 1. A single ion with 
mass M is trapped in a one-dimensional harmonic trap 
with the vibrational frequency v. \n) represents the n- 
th phonon state; fdand b are the phonon’s creation and 
annihilation operators. The internal states of the ion 
consist of three nondissipation states, | d), |r) and | g), and 
one dissipative state, |e), and the corresponding energies 
are Huiq, huj r , hoj g and hu > e , respectively. The lincwidth 
of state |e) is 2y. The transitions |e) -H- | g) and |e) O 
|r) are driven by two lasers with frequencies Wi, and 
wl 2 , wave vectors kL x and Vl 2 , and Rabi frequencies 9 r 


Hat = ~ Ag |e) (e| - A gr |r) (r\ 

+ (fig |e) (g| + fl r |e) (r| + fl M w Is) (d\ + H.c.). 

(7) 

represents the atomic internal DOF, and 

V = irjiflg |e) ( g\ (b + b *) + ir] 2 flr |e) (r| ( b + 6 t ) + H.c. 

( 8 ) 

describes the sideband transitions due to the two lasers. 
Here, the LD parameters are defined by 

vS: • ^ (9) 













If Qmw = A gr , we can rewrite the Hamiltonian of 
the atomic internal DOF 0 in a new representation 
{\D) ,\B) ,\+) ,\e)} 

H a t = — (A g ~ A gr ) |e) (e| + 2A gr |+) ( + | /i r,\ 

+ (Q+ |+) (e| + fls \B) (e| + H.c .), [W) 

where 

+ \/2 9 ’ (XI) 

1+) = ^(|$> + |d», 1 j 

= yUfTof (^9 l“) + l r )) ' 

It can be found that for the dark state \D), which is 


that is A + = A r — v, then EIT will arise between 
\D)\n) O |e) \n + 1} and |+}|n + l) O |e)|n+l), and 
the blue sideband heating is coherently canceled while 
only the cooling transition is left. As expected, FigE] 
shows a typical scattering spectrum of the cooling laser 
(the cooling transition |r) F4 |e)). The spectrum can 
be understood as consisting of two Fano-like structures, 
which is characteristic of EIT, and is controlled by the 
strong coupling laser field. Under the conditions EalTTl 
and EqEH all the heating transitions vanish, while the 
maximal red-sideband cooling is obtained. Therefore, the 
cooling cycle will start from the nth vibrational steady 
\D)\n), scatter over the excited state |e)|n — 1) via the 
red sideband transition, and dissipate back to the lower 
vibrational steady state \D)\n— 1), and the final double¬ 
dark state is 



where |—) = ^ (|g) — |0)), the state | D) decouples from 
the Hamiltonian (Hhd, which is the internal steady state, 
and the carrier transition vanishes due to the EIT effect. 
Then we can rewrite the sideband transition V with the 
representation (| D) , | B) , |+), |e}} as follows:, 

V = h]D^D\D) (e\(b + b^) + H.c. 

+ |+) (e| (b + fef) + H.c. (13) 

+ irjs^B \B) (e| (b + tf) + H.c. 


I'0) = H+ \D) |0) - ir] D n D |+) |1) + O (rj 2 ) , (18) 

which satisfies Hbd l^P) = 0. The final temperature can 
be evaluated by the mean phonon number of the double¬ 
dark state, yielding 


n ss — 


V 


D 


+ n 2 + 


otf). 


(19) 


According to (TTTT) and ED. 
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^ + 2 (^ + 1)0202 


+ O (v 4 ) ■ (20) 


with 


Vd = m - V2 ,Vb = 


Q' 2 g r)i + 2 C 1 2 t ]2 
202 + n 2 ■ 


O D = 




^202 + 02' 
(14) 

To understand the cooling dynamics, it is con¬ 
venient to transfer the Hamiltonian to a dressed 
representation!|+), |—), |r), \e)} , which is 


H at — A r |r)(r| + A_|—)(—| + A + |+)(+| 

+ ^(l+)(e| + |e)(+|) + 5|(|-)(e| + |e)H) 

+ O r (|r)(e| + |e)(r|). 


with 


A_ = A r 
A + = A r 


2A„ 


(15) 


(16) 


where A r , A+, A_ are corresponding detunings for tran¬ 
sitions 17') -H- |e), |+) - 0 - |e), |—) o |e), respective¬ 
ly. The transitions |—) O |e) and |r) -n- |e) have the 
same detunings to produce the dark state | D). Hence 
the carrier transition is eliminated, and the transitions 
\D)\n) O |e)|n + l) and |+)|ra + l) O |e) \n +1) re¬ 
main. If detuning A gr satisfies 


Agr ~ 2 ’ 


(17) 


If we choose 0 ff 0 r , the mean phonon number n ss will 
be much lower than the phonon recoil limit O (v 2 )- 



A f -A (in units of v) 


FIG. 2: (Color online) Absorption of cooling laser. As there 
are two EIT structures, the absorption vanishes at two values 
for the detuning A r of the cooling laser. Dashed lines mark 
the positions of these zeros and the absorption maximum. For 
appropriately chosen parameters, both the carrier transition 
and the blue-sideband heating are completely cancelled, while 
the red-sideband cooling is maximized. Here, the parameters 
are Q g = 7 / 7 / 2 , = y/20, A_ = 7 which represents the 

detuning for the transition |—) «-» \e). The other parameters 
are chosen to satisfy Ea. dlTI) and Ea. (l34l) . 

















III. ANALYTICAL TREATMENT 

The master equation is useful for analytically charac¬ 
terizing the cooling dynamics. Both the coherent dynam¬ 
ics and the dissipative nature of the excited level can be 
accounted for as follows: 

^p = £p = -i[H,p\+£ 1 p, (21) 

where the Liouvilian C y is for three dissipative channels, 
and 

C iP = J2 T* 

j=g,r,d 

~ |e)<e|p — p|e>(e|^, 

( 22 ) 

where 2 q j is the spontanous decay rate from |e) to 
\j) (j = 9 i r ,d), W (s) is the angular distribution of a 
spontaneous emission for a dipole transition, and kj = 
(w e - Wj) /c. 

In view of the small magnitude of 77 , the coupling be¬ 
tween the internal and external DOF (i.e., the laser cool¬ 
ing) is slow compared with the time scale of the internal 
atomic dynamics in the LDL. This situation allows adi¬ 
abatic elimination of the internal DOF, and thus we can 
obtain the rate equation of the motional state 

^Pn = {n+ l)A_p n+1 + nA + p n _i ^ 3 ) 
- [(n+l)A + + nA_\p n , 

where p n = (n\ p\n) is the population of the nth vibra¬ 
tional state and A± are the cooling and heating transition 
rates. 

Then the rate equation for the mean phonon number 
(n) is 

f t (n) = -(A.-A + )(n)+A + . (24) 

If A + < A_, the mean phonon number varies with 
time, as follows: 

(»} = «n) 0 -(n)Je-^ + (n) ss , (25) 

where ( n) 0 is the initial mean phonon number; (n) ss is 
the final mean phonon number; and W is the rate of the 
cooling process: 

(n)ss = a_- + A + +° (??2 )’ (26) 

W = A- — A + . (27) 

To obtain the cooling transition rates A±, we can de¬ 
fine r n ^„±i as the cooling and heat scattering ampli¬ 
tudes: 

F„^„±i =2 1 \{n±l\(D\W 1 G{E n )V\D)\n)\ 2 , (28) 


dsW(s)e l ^ X8 \j)(e\p\e)(j\e~ 


1-1 


where T n ^. n _i and T n ^,. n+ i describe the cooling and 
heating scattering process from the nth-vibration steady 
state \D)\n) to the (n ± l)th-vibration steady state 
\D)\n ± 1) via the sideband transition V. E n is the en¬ 
ergy of state \D)\n); the operator W y = |D)(e| describes 
the dissipation back to the internal dark state | D) from 
the excited state |e); and the function G(z) is defined by 


G(z) 


1 

z - H eS ’ 


(29) 


where the effective Hamiltonian H e g = — * 7 |e)(e|. 

Calculating EqQO we can get 


r„->„±i = 2^n 2 D r] 2 (n + 5±)v 2 


(Tv ~ 2A gr ) 2 

\f(Tv)\ 2 


(30) 


where <5+ = 1 ,S„ = 0, and 

f(x) = x 3 + (A g - A gr )x 2 

T (fl ^ 2 + I2 _|_ 2 + 2A 0 A ffr + 2A 2 gr )x + 2f2s“A gr 
+ i^(x — 2A gr )x. 

(31) 


Because = (n ± <5±)A±, the following holds: 

A± = (32) 


As expected, vanishes at A gr = — which means 
the carrier and the blue sideband heating are inhibit¬ 
ed,and the final mean phonon number n ss will be zero in 
the leading order of the expansion in the LD parameter 
77 . Then the corresponding cooling transition rate is 


A n — m—1 — 27D^)77 


Av 2 


(3 v 2 + 2vA g — 2 VL 2 b — 0+) + A r ) 2 i' 2 

( 33 ) 

We can obtain the maximal cooling rate at 


3 v 2 + 2vA g - 2n 2 B -n 2 + =o. 


(34) 


The maximum W is 


W = 


2£l 2 D p 2 

7 


(35) 


IV. NUMERICAL RESULTS AND DISCUSSION 

We have described the model and obtained the analyt¬ 
ical result of the cooling process. In this section, we use 
the quantum toolbox [29| for numerical simulation of the 
master equation Eq[3] A m Yb + ion is regarded as an 
excellent candidate for lager-scale quantum computation 
with qubits encoded in two hyperfine states H3, [30- 
[34j |. Our scheme can be experimentally implemented 
in the m Yb + ion system, where the energy levels of 
|r), 1 3 ), \d) and |e) can be realized by choosing | 2 5 i/ 2 , F = 









FIG. 3: (Color online) Numerical simulation of cooling a 
1 , 1 Yb + . The parameters are r] g = —r} T = 0.05, 2y = 20u, 
7 g ~ 7 r — 7 d ~ 10i//3, fig — lOu, fl r — If, H.MW — Ag r 
5 Ag = 74.5^. The analytical predictions are shown as 
dashed lines.The solid lines correspond to the numerical re¬ 
sults. The final phonon number is much less than the recoil 
limit 0( 10 -2 ). 


1 ,m F = —1), \ 2 Si/ 2 ,F = 1 ,m F = 0),\ 2 S 1 /2,F = 
1 ,mj 7 = 1 ) and \ 2 P\/ 2 ,F = 0 , TOf = 0 ), respectively.We 
can use two 369-nm laser beams to couple the transition 
\ 2 S l/2 ,F = 1 ,m F = -1) o \ 2 Pi/ 2 ,F = 0 ,m F = 0) 
and \ 2 S 1 / 2 ,F = 1 ,m F = 0) \ 2 Pi/ 2 ,P = 0 ,m F = 0), 

and use an acousto-optical modulator (AOM) driver to 
tune the detuning difference between these two transi¬ 
tions. To drive the transition \ 2 Si/ 2 ,F = 1 ,mp = 0) O 
| 2 5 i/ 2 > P = 1 ) m F = 1 ), a magnetic field can be used to s- 
plit these two generate levels, and then the transition can 
be driven with a microwave field.The branching ratios out 
of the excitation state \ 2 Pi/ 2 ,F = 0 , tof = 0 ) to the three 
hyperfine states. \ 2 Si/ 2 ,F = l,mj? = — 1), | 2 S' 1 / 2 , F = 
1 ,tof = 0) and \ 2 Si/ 2 ,F = 1 ,tof = 1) are the same, 
that is 7 g = 7 r = 7 d- 

As shown in Fig. [3l the analytical prediction of the 
cooling dynamics in Eq[25] is consistent with the nu¬ 
merical simulation, and the final phonon number of our 
scheme can reach the order of 0( 10 -4 ). As t]d = 0.1, the 
corresponding recoil energy is approximately ~ 

O( 10 ~ 2 ), and then the final temperature in our scheme 
can be much lower than the recoil energy for Q g fl r . 

The dependence of the final cooling result on the Rabi 
frequencies fl g , f2 r is shown in Fig. [Hand Fig. [5l respec¬ 
tively. In Fig. 01 with the fixed f2 r and other parameters 
chosen under optimal conditions m and (l34l) . the final 
mean phonon number decreases with increasing Sl g . In 
contrast, as the Rabi frequency is fixed, the final cool¬ 
ing result worsens with increasing which is shown in 
Fig. [5] In the two figures, the numerical simulation a- 
grees with the analytical prediction Ea. (l2Ul) . The final 
phonon number increases by increasing the strength of 
f l g while it decreases with increasing f l r because they 
affect the final double-dark according to Eq. [H]. 

Our cooling scheme is robust with respect to the fluc- 



FIG. 4: (Color online) Numerical simulation of the final 
phonon numbers as a function of fl 9 . The parameters are 
r/ g = —7] r = 0.05, 2y = 20^, 7g = 7r- = 7d = 10^/3, fT = lv, 
Omw = A gr = —0.5 f, A g are chosen under the optimal con¬ 
dition Eq. m- The dashed line is obtained by analytical pre¬ 
diction while the solid line corresponds to numerical results. 
The cooling result gets better with increasing fig. 



FIG. 5: (Color online) Numerical simulation of the final 
phonon numbers as a function of Q r - The parameters are 
r/g = —Vr = 0.05, 27 = 20f, 7g = y r = 7d = 10i//3, = IO 17 

O mw = A gr = —0.5^, A g are chosen under the optimal 
condition Eg. (l34II .The dashed line is obtained with analyti¬ 
cal prediction while the solid line corresponds to numerical 
results. The curves show the final phonon number decreases 
with increasing f l r . 


tuation of the Rabi frequencies. As depicted in Fig. [ 6 ] or 
Fig. [7] when the fluctuation of the Rabi frequency FImw 
or f l g is approximately 20 %, the final mean phonon num¬ 
ber varies by the order of 10 -3 . Therefore, the cooling 
scheme is insensitive to the fluctuation of the microwave 
power and the laser optical intensity. 

Our scheme is more suitable for the ions with hyper¬ 
fine structures such as m Yb + and m Cd + . For pre¬ 
vious methods based on the three-level A-configuration, 
in experimental implementation with special ions such 
as 40 Ca + and 24 Mg + , unwanted upper levels result in 
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FIG. 6: (Color online) Numerical simulation of the final 
phonon numbers as a function of Qmw- The parameters are 
Tj g = —r/r = 0.05, 27 = 20u, 7 g = -fr = 7 d = l(fo/3, fig = IO 17 
fl r = 1 v, A gr = —0.5 v, A g = 74.5 v. The corresponding 
optimal value of flMW from the analytical prediction is at 
flMwo = — 0 . 5 za The cooling result is robust against the 
fluctuation of the microwave driver’s power flMW- 


x10' 4 



FIG. 7: (Color online) Numerical simulation of the final 
phonon numbers as a function of fl g . The parameters are 
Vg = —Vr = 0.05, 27 = 20v, 7 9 = 7 r = 7 d = 10zx/3, fl r = 1 v, 
flMW = A 9 r = —0.5ix, A g = 74.5^. The corresponding opti¬ 
mal value of fig from the analytical prediction is at fl g = 10 v. 
The cooling result is robust against the fluctuation of fl g that 
depends on the optical laser intensity. 

decreasing the cooling performancefit! [l9j. whereas for 
1 (1 Yb + and other parameters here, the unwanted hyper- 
fine structure \ 2 Pi/2 ,F = 1) which includes three hyper- 
fine excitation states in radiant 2y = 19.7 MHz has a nat¬ 
ural 2.1 GHz gap with the level \ 2 Pi/ 2 ,F = 0, tof = 0), 
and only \ 2 P 1 / 2 ,F = 1 ,m F = -1), \ 2 Pi/ 2 ,F = 0 ,m F = 

1) can lead to a weaker scattering rate R ~ | x (%-) 2 x 
2y = 0.32 KHz compared to the cooling rate W = 2 KHz 
according to Eq. [H]. Hence the cooling scheme here is ef¬ 
ficient. Moreover, we only require one laser to implement 
Doppler cooling and our cooling scheme simultaneously 
instead of Raman sideband cooling. 


Recently, large-scale quantum-computer architecture 
are more and more attractive [24h281 130)-(40| , and scaling 
can potentially be achieved by storing ions in multizone 
arrays where information is moved in the processor by 
physically transporting the ions ,3jl, 36] or teleporting 
[371 ] . It is important to recool the qubit ion because ion- 
s can be heated by ambient noisy electric fields and/or 
during ion transport j3o|. Experimentally, “refrigerant” 
ions can be used to cool the qubit ions sympathetical¬ 
ly 0 ES HI; HU HM3 ■ Microwave can be applied via 
waveguides (structures that can guide radiation) that are 
part of the chip on which the ion trap is integrated [27| , 
and thus does not require alignment, while laser beams 
must be carefully aligned to interact with the trapped 
ions [26}. Moreover, it is easier to generate and control 
microwave fields which are more robust to amplitude, 
phase, frequency and polarization noise compared to the 
use of lasers [2J, |25| , and these parameters are difficult 
to control in integrated optics |28j. Hence, compared to 
double-EIT cooling, our cooling scheme uses a microwave 
driver instead of one laser beam to ease the experimental 
condition for cooling “refrigerant” ions in multizone ar¬ 
rays, and make it more feasible to realize the double-dark 
state cooling on the ion chip firstly. 


V. CONCLUSION 

In conclusion, we have shown a new efficient double¬ 
dark state cooling scheme for a trapped ion system us¬ 
ing the EIT effect combined with microwave Stark shift. 
Microwave coupling is used to dress the ion levels and 
obtain two EIT structures, providing a double dark s- 
tate. In leading order of the Lamb-Dicke expansion, this 
suppresses all heating excitations. As a consequence, the 
final temperature in this scheme can be much lower than 
the recoil energy. It is robust against fluctuations of the 
microwave power and the laser intensities. Compared 
with other dark-state cooling methods, this scheme is 
better for the experimental realization of the cooling of 
some specific trapped ions. 


ACKNOWLEDGMENTS 

We thank Y. Lin and Q. Yin for helpful comments 
on the manuscript. This work was supported by the 
National Natural Science Foundation of China (Nos. 
61327901, 61490711, 11274289, 11325419, 11474268, 
11304366, 11404319, 91421111, 61205108), the National 
Basic Research Program of China (No. 2011CB921200), 
the Strategic Priority Research Program (B) of the Chi¬ 
nese Academy of Sciences (Grant No. XDB01030300), 
the Fundamental Research Funds for the Central Uni¬ 
versities (No. WK2470000011, No. WK2470000018), 
the China Postdoctoral Science Foundation (Grant Nos. 
2013M531771 and 2014T70760), Foundation of Science 
and Technology on Information Assurance Laboratory 




(No. KJ-14-001); 


[1] Cirac, J. I. and P. Zoller, Phys. Rev. Lett. 74, 4091 
(1995). 

[2] C. Monroe, D. M. Meekhof, B. E. King, W. M. Itano, 
and D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995) 

[3] Childs, A. M. and I. L. Chuang, Phys. Rev. A. 63, 012306 

( 2000 ). 

[4] P. O. Schmidt, T. Rosenband, C. Langer, W. M. Itano, J. 

C. Bergquist, and D. J. Wineland, Science. 309, 749-752 
(2005) 

[5] F. Diedrich, J. C. Bergquist, W. M. Itano, and D. J. 
Wineland, Phys. Rev. Lett. 62, 403 (1989) 

[6] C. Monroe, D. M. Meekhof, B. E. King, S. R. Jefferts, 
W. M. Itano, D. J. Wineland, and P. Gould, Phys. Rev. 
Lett. 75, 4011 (1995) 

[7] C. Roos, T. Zeiger, H. Rohde, H. C. Nagerl, J. Eschner, 

D. Leibfried, F. Schmidt-Kaler, and R. Blatt, Phys. Rev. 
Lett. 83, 4713 (1999) 

[81 Stefano Zippilli and Giovanna Morigi, Phys. Rev. A. 72, 
053408 (2005) 

[9] Marc Bienert and Giovanna Morigi, New. J. Phys. 14, 
023002 (2012) 

[10] Shuo Zhang, Qian-Heng Duan, Chu Guo, Chun-Wang 
Wu, Wei Wu, and Ping-Xing Chen, Phys. Rev. A 89, 
013402 (2014) 

[11] Giovanna Morigi, Jurgen Eschner, and Christoph H. Kei¬ 
tel, Phys. Rev. Lett. 85, 4458 (2000) 

[12] J. Evers and C. H. Keitel, Europhys. Lett. 68, 370 (2004) 

[13] A Retzker and M B Plenio, New. J. Pliys. 9, 279 (2007) 

[14] J. Cerrillo, A. Retzker, and M. B. Plenio, Phys. Rev. 
Lett. 104, 043003 (2010) 

[15] Giovanna Morigi, Phys. Rev. A. 67, 033402 (2003) 

[16] C. F. Roos, D. Leibfried, A. Mundt, F. Schmidt-Kaler, J. 
Eschner, and R. Blatt , Phys. Rev. Lett. 85, 5547 (2000) 

[17] Shuo Zhang, Chun-Wang Wu, and Ping-Xing Chen, 
Phys. Rev. A. 85, 053420 (2012) 

[18] Shuo Zhang, Jian-Qi Zhang, Qian-Heng Duan, Chu Guo, 
Chun-Wang Wu, Wei Wu, and Ping-Xing Chen, Phys. 
Rev. A. 90, 043409 (2014) 

[19] Y. Lin, J. P. Gaebler, T. R. Tan, R. Bowler, J. D. Jost, 
D. Leibfried, and D. J. Wineland, Phys. Rev. Lett. 110, 
153002 (2013) 

[20] Keyu Xia and Jorg Evers, Phys. Rev. Lett. 103, 227203 
(2009) 

[21] J. Q. Zhang, S Zhang, J. H. Zou, L. Chen, W. Yang, Y. 
Li, and M. Feng, Optics Express 21, 29695-29710 (2013) 


[22] Michael Fleischhauer, Atac Imamoglu, and Jonathan P. 
Marangos , Rev. Mod. Phys. 77, 633 (2005) 

[23] A. Albrecht, A. Retzker, C. Wunderlich and M. B. Plenio, 
New. J. Phys. 13, 033009 (2011) 

[24] C. Ospelkaus, U. Warring, Y. Colombe, K. R. Brown, J. 
M. Amini, D. Leibfried, D. J. Wineland, Nature. 476, 
181C184 (2011) 

[25] D. P. L. Aude Craik, N. M. Linke, T. P. Harty, C. J. 
Ballance, D. M. Lucas, A. M. Steane, D. T. C. Allcock, 
Applied Physics B. 114, 3-10 (2014) 

[26] W. K. Hensinger, Nature. 476, 155 (2011) 

[27] M. D. Hughes, B. Lekitsch, J. A. Broersma and W. K. 
Hensinger, Contemp. Phys. 52 505C29 (2011) 

[28] G. R. Brady, A. R. Ellis, D. L. Moehring, D. Stick, C. 
Highstrete, K. M. Fortier, M. G. Blain, R. A. Haltli, A. A. 
Cruz-Cabrera, R. D. Briggs, J. R. Wendt, T. R. Carter, 
S. Samora, S. A. Kemme, Appl. Pliys. B Lasers Opt. 
103(4), 801 (2011). 

[29] Sze.M.Tan, J.Opt. B. 4, 424 (1999) 

[30] C. Monroe and J. Kim, Science 339, 1164 (2013). 

[31] C. Monroe, R. Raussendorf, A. Ruthven, K. R. Brown, 
P. Maunz, L.-M. Duan, J. Kim, Phys. Rev. A 89, 022317 
(2014). 

[32] T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. 
Monroe, and J. L. OBrien, Nature 464, 45 (2010). 

[33] R. C. Sterling, H. Rattanasonti, S. Weidt, K. Lake, P. S- 
rinivasan, S. C. Webster, M. Kraft ,and W. K. Hensinger, 
Nat.Commun.5, 3637, (2013) 

[34] Ch. Piltz, Th. Sriarunothai, A. F. Varn, and Ch. Wun¬ 
derlich, Nat.Commun.5, 4679 (2014). 

[35] D. J. Wineland, C. Monroe, W. M. Itano, D. Leibfried, 
B. E. King, and D. M. Meekhof, J. Res. Natl. Inst. Stand. 
Technol. 103, 259 (1998). 

[36] D. Kielpinski, C. Monroe, and D. J. Wineland, Nature 
(London) 417, 709 (2002). 

[37] D. Gottesman and I. L. Chuang, Nature (London) 402, 
390 (1999). 

[38] J. D. Jost, J. P. Home, J. M. Amini, D. Hanneke, R. 
Ozeri, C. Langer, J. J. Bollinger, D. Leibfried, and D. J. 
Wineland, Nature (London) 459, 683 (2009). 

[39] J. P. Home, D. Hanneke, J. D. Jost, J. M. Amini, D. 
Leibfried, and D. J. Wineland, Science 325, 1227 (2009). 

[40] D. Hanneke, J. P. Home, J. D. Jost, J. M. Amini, D. 
Leibfried, and D. J. Wineland, Nat. Phys.6, 13 (2009). 



